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I 

Abstract 

We study five-dimensional black holes obtained by compactifying M theory on Calabi- 
Yau threefolds. Recent progress in solving topological string theory on compact, one- 
parameter models allows us to test numerically various conjectures about these black 
holes. We give convincing evidence that a microscopic description based on Gopakumar- 
Vafa invariants accounts correctly for their macroscopic entropy, and we check that highly 
nontrivial cancellations -which seem necessary to resolve the so-called entropy enigma in 
the OSV conjecture- do in fact occur. We also study analytically small 5d black holes 
obtained by wrapping M2 branes in the fiber of K3 fibrations. By using heterotic/type II 
duality we obtain exact formulae for the microscopic degeneracies in various geometries, 
and we compute their asymptotic expansion for large charges. 
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1 Introduction 

String theory can provide in many situations a precise microscopic description of su- 
persymmetric black holes which reproduces for large charges the Bekenstein-Hawking 
entropy [39J. Degeneracies of microstates that are highly protected by supersymmetry 
are often counted by mathematically well understood quasitopological quantities related 
to the compactification manifold. For example, the computation of microstates of the 
D1-D5 system is encoded in the elliptic genus of a symmetric product of a hyperKahler 
manifold (see [HI E] for a review of these computations). 

A very challenging class of black holes in string theory is obtained by compactifying 
M theory on a Calabi-Yau manifold X with generic SU(3) holonomy. These are five 
dimensional black holes, which are characterized by a membrane charge Q G ^(A, Z) 
and an angular momentum m. It was proposed in [28] that the microscopic entropy of 
these black holes is accounted for by BPS states associated to M2 branes wrapping the 
cycle Q and with left spin m = 2ji in five dimensions. Accoording to the proposal of [28] 
their degeneracies are encoded in the five dimensional supersymmetric index 

I(a, (3) = Tr(-l) 2jL+2jR exp(-aj L - (3H). (1.1) 
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The information captured by this index can be extracted from the all-genus expansion of 
the holomorphic free energy of the topological string, computed at the large radius point 
of .V EH- 

oo 

limF(t,t,g s ) = J2g 2 r 2 F 9 (t)- (1-2) 

t— >oo ^ ' 

3=0 

In this identification the BPS degeneracies are mapped to the Gromov-Witten invariants 
of genus g holomorphic maps. Since the topological string could not be solved on com- 
pact Calabi-Yau threefolds at higher genus, progress in understanding the microscopic 
degeneracies in the general case was very limited. On the boundary of the Kahler cone 
the problem might reduce effectively to a counting problem on a two complex dimensional 
surface, which is mathematically simpler, but the situation is also physically more degen- 
erate. When the compactification manifold is X = X2 x T 2 , where X2 = T 4 or K3, one 
obtains the five dimensional black hole solutions constructed in [5], and the microscopic 
degeneracies are encoded in the elliptic genus of symmetric products of X 2 . 

In this paper we study the microscopic counting proposed in [2H] in two different situ- 
ations, by using numerical and analytic methods. First of all, we consider 5d black holes 
obtained by compactifying M-theory on the one-parameter Calabi-Yau spaces studied in 
[26] . This explores a generic direction in the Kahler cone and allows to describe generic 
5d black holes, which have non-vanishing classical horizon area and can carry spin. In 
[2H] significant progress was made in solving the topological string on compact Calabi- 
Yau threefolds. By combining the holomorphic anomaly of [I] with modularity properties 
of the topological string partition function Z = exp(F), effective action arguments, and 
Castelnouvo theory, it was possible to calculate the topological string free energy up to 
genus 53. 

In order to make contact with black hole physics on the (super)gravity side, one 
has to obtain the asymptotic expansion of the microscopic degeneracies for large charge 
Q and Q ^> m. For fixed g the expansion of F g (t) around large radius is convergent 
and under analytic control by mirror symmetry. In contrast, the genus expansion in 
fll.2p is expected to be asymptotic, as in noncritical string theories [38] (see [33] for a 
recent discussion of this issue in the context of topological strings). To obtain a large 
Q expansion for the degeneracies of the (Q, m) states one needs information at genus 
g ~ Q 2 and is hence facing the issues of the behavior of string theory at large genus. 
Although we don't have enough control of the degeneracies to obtain analytical results 
on the large charge expansion, the situation is suited to a numerical study by using 
the method of Richardson transforms^. This method merely relies on the knowledge of 
the form of the series and makes it possible to extract its coefficients from the value 
of the degeneracies at finitely many points. The analysis is complicated by the fact 
that the large charge expansion of the degeneracies is an asymptotic expansion, but we 
find that the Richardson transforms converge rapidly to the expected macroscopic values 

^-For sub-sub leading terms we use the Pade approximation. 
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for the asymptotic coefficients. To estimate its accuracy we sample over the thirteen 
Calabi-Yau, which have a sufficently wide variety of topological data. Using this sample 
we can conclude that, given our present data of the higher genus instanton expansion, 
the leading coefficient of the asymptotic expansion is correct within 2% and the first 
subleading one within 12%. With this information at hand, we give convincing evidence 
that the topological string accounts correctly for the entropy of 5d spinning black holes, 
as conjectured in |28^l . 

Some aspects of the genus expansion (11.21) are much better understood in terms of 
D-brane invariants like Gopakumar-Vafa or Donaldson-Thomas invariants, rather than 
Gromov— Witten invariants. In particular, for a given charge Q G H 2 (X, Z) and Q ^ one 
gets the complete genus information from a finite number of GV invariants. We use the 
results for F g in [26] to obtain precise information on the Donaldson-Thomas invariants 
of the one-parameter models. This allows us to study numerically the scaling exponent 
k considered in [12J (and defined below in (14 .7p and (I4.8P ). which governs the growth 
of the Donaldson-Thomas invariants under rescalings of the charges. Our numerical 
study indicates that k = 2. As argued in [T2], this value indicates that highly nontrivial 
cancellations occur between the contributions to the Donaldson-Thomas invariants, which 
in turn seem necessary to resolve the so called entropy enigma [12] in the OSV conjecture 

The second class of black holes we study has a different flavor. These are 5d black holes 
which are obtained when the Calabi-Yau is a K3 fibration and the charge Q is restricted 
to the K3 fiber. Their classical horizon area is zero (i.e. they are small black holes) 
and have no spin. By using heterotic/type II duality one can obtain analytic formulae 
for the F g amplitudes at all g [2j [3U [3QI E2l [22], and from them one can extract the 
exact microscopic degeneracies for the corresponding small 5d black holes. Of course, as 
explained for example in [T2], the most delicate aspects of 5d spinning black holes, as 
well as of the OSV conjecture, cannot be tested with small black holes. This reflects the 
fact that the Gromov- Witten theory of K3 fibers (which is closely related to the theory 
of Hilbert schemes) is much simpler than the Gromov- Witten theory of generic Calabi- 
Yau manifolds. However, having an exact microscopic counting might be important in 
understanding some detailed aspects of the entropy. As in the 4d case considered in [5] , the 
5d degeneracies are closely related to modular forms, but one cannot use the Rademacher 
expansion featured in [T5], [9]. We find however an exact asymptotic expansion for the 
microscopic degeneracies in powers of the inverse charge (albeit corrected by terms which 
are exponentially suppressed for large charges). The leading term of the asymptotics is 
in agreement with a macroscopic computation using the 4d/5d connection of [20] and the 
4d attractor equations for a D6/D2 system. 

The organization of this paper is as follows. In section 2 we review the macroscopic 
and microscopic computation of the entropy for 5d spinning black holes. In section 3 

2 For a recent study of this question by using an approach totally different from ours, see [24] . 
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we analyze numerically the asymptotic properties of the degeneracies for the Calabi-Yau 
manifolds studied in [26]. In section 4 we study the asymptotic properties of Donaldson- 
Thomas invariants to address the entropy enigma of [12]. In section 5 we study small 
black holes in K3 fibrations and compute their degeneracies as well as the asymptotic 
expansion. Finally, in section 6 we list some conclusions and open problems. 



2 Microscopic and macroscopic counting for 5d black 
holes 

2.1 Macroscopic description 

Let us start with the macroscopic description of black hole entropy. We will consider 
5d black holes obtained by compactifying M theory on a Calabi-Yau threefold X, and 
characterized by a charge Q e H 2 (X,Z,) and SU{2) L C 5*0(4) angular momentum to. 
We will introduce a basis S A for H 2 (X, Z), where A = 1, • • ■ ,b 2 (X), as well as a dual 
basis tja for H 2 (X). With respect to the E" 4 basis, the charge Q will be given by a set of 
integers Qa- The classical entropy of these black holes, denoted as So, is one quarter of 
the horizon area 



S = 27r v / Q 3 -m 2 , (2.1) 

where Q is the graviphoton charge of the black hole. This charge is related to the mem- 
brane charge Q by the attractor mechanism in five dimensions |19| . which states that 



Q 3/2 = \c ABC y A y B y c , (2.2) 
o 

where 

\c ABC y B y c = Qa, (2.3) 

and 

Cabc = / VaAVbAVc (2.4) 
Jx 

are the triple intersection numbers of X. For one-parameter models, the membrane charge 
will be identified with the degree d of the holomorphic map in topological string compu- 
tations, and we will denote the single intersection number by Cabc — K- From the above 
equations it follows that 

Q=l7T-M- (2-5) 



9k / 

There is a correction to the black hole entropy from the R 2 term of the supergravity 
effective action, which we denote as Si for convenience. The R 2 term correction to the 
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black hole entropy scales like in the large charge limit, and was computed in [23] by 
using Wald's formula [JT] for the R 2 in 5d. The result reads, 



(2.6) 



where 

Y A = ^ (2-7) 

and 

c 2 a= [ c 2 (X)A VA . (2.8) 
Jx 

For m = this formula has been rederived in [SI [I] by using the full 5d SUGRA action. 
In the one-parameter case, this correction reads 

Besides the corrections that we have considered, there are well known correction terms 
in the low energy effective action of the form F 2a ~ 2 R 2 , g > 2, where F is the graviphoton 
field strength. The leading contribution comes from a classical term, which is the contri- 
bution of the constant map from a genus g world-sheet to the Calabi-Yau manifold. It is 
of the form, 

d gX (2.10) 
where \ is the Euler number of the Calabi-Yau 3-fold and 

da Ag{2g-2){2g-2)V ^ 

We denote the correction to black hole entropy due to the F 2g ~ 2 R 2 term as S g . We can 
roughly estimate the correction for non-spinning black holes m = as follows. 

The graviphoton charge is the integral of its field strength over the horizon of black 
hole, 

Q~ / F (2.12) 

J Horizon 

Since the area of the horizon scales like A ~ Qa, the graviphoton field strength goes like 

F ~ Q-h (2.13) 
The R 2 term contributes a factor Q _1 in Wald's formula for the black hole entropy, and 

3 

taking into account also the factor of horizon area A ~ Qa ; we find the scaling behavior 
of the F 2g ~ 2 R 2 term correction to black hole entropy to be 

S g ~ xQ 1 ' 9 (2.14) 
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where we have included the Euler number from (12.101) . The constant of proportionality in 
(12.141) is now universal, and independent of specific Calabi-Yau geometries and the black 
hole charge. We will be able to make a rough test of (12.141) for the genus 2 case, which is 
the sub- sub- leading correction in the large Q limit. 

There are other worldsheet instanton corrections to the low energy effective action 
that can be computed also by topological strings. However, these terms are exponentially 
small in large charge Q limit where the supergravity description is valid, and are much 
suppressed compared to the Q~ l power corrections in (I2.14p . In this paper we will not 
need to consider these world-sheet instanton corrections in the macroscopic description of 
the black hole entropy. Interestingly these world-sheet instanton corrections are closely 
related to the BPS states that we will count in the microscopic description of the black 
hole entropy. 



2.2 Microscopic description 

Microscopically, a 5d black hole with membrane charge Q e H 2 (X,'Z) is engineered by 
wrapping M2 branes around the two-cycle Q. This leads to a supersymmetric spectrum 
of BPS states in 5d which are labeled by Q and by their spin content (Jl, ]r)- As argued 
in [21], in order to obtain an index one has to trace over jn (with an insertion of (— l) 2jii ). 
The resulting spectrum for a membrane charge Q can be represented as 

9 
r=0 

where 




(2.16) 



encodes the spin content ji, and Uq are the Gopakumar-Vafa invariants [21 J. Notice that 
in (12.151) . the sum over r is finite and the highest spin g appearing in the sum depends 
on the membrane charge Q. In other words, for any given Q there are only finitely many 
g for which the rig are nonzero. 

We can now write down a generating function for the supersymmetric degeneracies of 
BPS states with membrane charge Q, keeping track of their left spin ji, by computing 

X)n(Q,m) = J>* Q (-1) 2 V-. (2.17) 

m Q 

Using the decomposition (I2.15P one finds 
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where m = 2ji- In [28J it was proposed that this quantity gives the microscopic degen- 
eracies for a spinning 5d black hole of charge Q and spin J = m. The computation of 
these degeneracies reduces then to the computation of the Gopakumar-Vafa invariants 
rig. The most effective way to determine these is by computing the genus r amplitudes 
F r of topological string theory on X. As shown in [2T], the total free energy 

oo 

F(t,g s ) = J £ F r( t )9s r - 2 (2-19) 

r=0 

can be written in terms of the Gopakumar-Vafa invariants as 

oo oo ^ , U \ 2r— 2 

n*,*) = E E E^K 2sin 2 ) e " fcQ '- (2 - 20) 

r=0 QeH 2 (X) k=l ^ ' 

This means, in particular, that one can obtain the riq for all Q by knowing Fq, ■ ■ ■ ,F r . 
The black hole entropy is given by logarithm of the number of microstates 

S{Q,m)=log{tt{Q,m)). (2.21) 

This should agree with the macroscopic result in the large charge limit Q ^> 1 and Q ^> m. 

As explained in [28J , this proposal for the microscopic counting of states of 5d black 
holes can be regarded as a generalization of the elliptic genus, which computes BPS 
degeneracies of the D1-D5 system. Indeed, if one considers M theory compactified on 
X = K3 x T 2 , the generic M2 brane charge in this compactification is 

Q=[C]+n[T\ neZ, (2.22) 

and C is a curve in K3. By standard dualities this system can be related to type IIB 
on K3xS x with D-brane charge [C] and M units of momentum around S 1 , which is a 
close cousin of the D1-D5 system. As it is well known (see for example [14]). the BPS 
degeneracies of this system can be computed from the elliptic genus of the symmetric 
product of K3. Let 

x(K3;g,y)=EcK^)?V (2-23) 

m,i 

be the elliptic genus of K3. The generating function of elliptic genera of the symmetric 
product S k K3 

X (S P K3; g,y) = E x(S k K3; q, y)p k = £ c(k, n, m)p k q n y m (2.24) 

k=0 k,n,m 

can be computed from the coefficients in f!2.23j) in terms of an infinite product [16J 

X (S p K3;q,y)= ]J (l-p^V)" W) . (2.25) 

n,m>o/ 
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56 


8 


-^3,2,2 (1 ) 
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-128 


64 


16 
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-156 


48 


6 


Y 4 , 4 (l 4 ,2 2 ) 


-144 


40 


4 


X 6 , 2 (l 5 ,3) 


-256 


52 


4 


Y M (1 3 ,2 2 ,3) 


-156 


32 


2 


X 6 , 6 (l 2 ,2 2 ,3 2 ) 


-120 


32 


1 











Table 3.1: The sample of 13 one-parameter complete intersection CYs in weighted pro- 
jective spaces. A CICY of degree d\, . . . , c4 in weighted projective space F l ~ 1 (wi, . . . , w{) 
is denoted X dl> ...^(wi, . . . , wi), i.e. weights w with repetition m are abbreviated by w m . 
X = f x c 3 * s the Euler number, k is the triple intersection number, and c 2 • t] is defined in 

The supersymmetric degeneracies of BPS states for this system are then given by the 
coefficients of the expansion in (12.241) , 

Q(Q, m) = c(Uj 2 + 1, n, rnj , (2.26) 
where Q is of the form (I2.22p . One can show that, for large charges [T7J [14"] . 



log n(Q,m) ~ 2vr^-C 2 - m 2 . (2.27) 

It is easy to check that this is precisely the macroscopic entropy (12. ip computed for 
K3x T 2 . Of course, the degeneracies (12.181) are in general much more difficult to compute, 
since they correspond to black holes with only M = 1 supersymmetry in 5d. 

3 One— parameter models 

3.1 Topological strings on one— parameter models 

In [26] the topological B model was integrated on thirteen one-parameter Calabi-Yau 
spaces which can be realized as hypersurfaces or complete intersections in (weighted) 
projective spaces. We have listed these spaces and some of their topological data in table 
13.11 These data are the intersection numbers Cabc — K > the second Chern classes c 2 , and 
the Euler numbers \. They are needed for computations of the macroscopic entropy. 

The complex moduli space of these threefolds is M. = P 1 \ {oo,l,0}, and the three 
special points are the large radius degeneration point, a conifold point and a further point 
either of finite (Gepner point) or infinite branching order. The modular group Tx £ 
SP(4, Z) can hence be generated e.g. by the large radius and the conifold monodromies. 



S 



The conceptual obstacle in integrating the B-model holomorphic anomaly [I] is the 
holomorphic ambiguity which arises in each integration step. Invariance of the topological 
string amplitudes under Fx and effective action arguments, which govern the behaviour of 
the genus g amplitudes at special points, restrict the ambiguity to 3g — 3 unknowns |28j . 
By using a refined effective action analysis, which gives rise to the "gap condition" at 
the conifold, regularity at the orbifold, and Castelnouvo's bound for the Gromow-Witten 
invariants at large radius, it is possible to fix the unknows, and one can calculate the free 
energy of the topological string to arbitrary degree and up to genus 12 — 53. 

Instead of using the generic solution of holomorphic anomaly equation suggested 
by the world-sheet degenerations [lj we use the constraints of T x on the topological 
string amplitudes directly when integrating the holomorphic anomaly equations genus by 
genus [42] [22] . This results in an alogarithm, which constructs the genus g amplitudes 
as weight 3g — 3 polynomials over a ring of three an-holomorphic- and one holomorphic 
modular objects of weight (1,2,3, 1). As a consequence the number of terms in the F g 
grows polynomial with g and not exponentially as in the approach of [4^|. 

The approach of views the topological string partition function as a wave func- 
tion over H 3 (X, R). Choices of polarization are necessary in order to expand the effective 
action at different points in the moduli space Ai, in appropriate local holomorphic coordi- 
nates. Most of the black hole issues that we will discuss involve the degeneracies extracted 
from the topological string at the large radius limit. Therefore we will discard for now 
most of the global information and focus only the holomorphic limit of the topological 
string partition function at this limit, where it encodes the degeneracy of bound states of 
a single D6 brane and arbitrary D2-D0 branes. 



3.2 Static black holes 

We first consider the case of non-spinning black hole J = m = and denote = Q(d, 0). 
The entropy formula including the first few orders f!2.ip . (12. 9p . (12.141) is in this case 



S = b d* + b 1 d* + -J r + O — , (3.1) 

G?2 \d2 J 

where the first two coefficients are 

b = —=, 6i = 3.2 
3^2^ 4^2^ 

and we have used the graviphoton charge relation (12. 5p . The coefficient b\ of the sub- 
leading term is consistent with the results in [4"0l M, d] , where it was observed that the b\ 
can be produced by a shift of the 2-brane charge 

d^d+— (3.3) 



Nevertheless the limiting factor in advancing to higher g is presently not the ambiguity but the com- 
puting time. The reason is that the numerators in the coefficents of the polynomials grow exponentially. 
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in the leading term. 

To compare with the microscopic counting we define the following function 



m = (3.4) 

The macroscopic black hole entropy predicts that the large order behavior of f(d) is 

/(rf) = 6o + 7 + l + "" (3 ' 5) 

Since we have available only the values of f(d) for positive integer d up to a finite degree, it 
is appropriate to use well-known numerical methods to extrapolate the asymptotic value 
b . From the form of the sub- leading corrections in (I3.5p . it is appropriate to use the 
Richardson extrapolation method (see for example [3]). 

The basic idea of this numerical method is simple. To cancel the sub-leading correc- 
tions in (13.51) up to order l/d N , one defines 

One can show that if the perturbation series (13.51) truncates at order l/d N , the expression 
(13. 6p will give exactly the asymptotic value bo- Ideally ,the larger N and d are, the closer 
A(d, N) is to the asymptotic value. But due to the limitation of our data, the sum d+ N 
must not exceed the maximal degree d max of the topological string computations. 

Fig. [T] shows the convergence of the leading terms in f(d) and of the Richardson 
transforms A(d,N), N = 2,3,4 for the quintic and the bi-cubic. It is obvious from the 
two examples in Fig. [T]that the Richardson method improves impressively the convergence 
of the series, i.e. it provides a model independent and consistent scheme to supress the 
subleading corrections. Using iV = 2 — 4 is good enough for our purpose of estimating 
the asymptotic value. We conduct the analysis for all 13 models using N = 3 and the 
maximal degree available. The results are summarized in Table I3.2[ and are in very good 
agreements with the expected asymptotic values bo in (I3.2p . More detailed results on all 
the analysis carried out in this paper can be found in a script and in a data base at |43j . 

We can further extract the sub-leading coefficient b\ from the data. Define 

hid) = (f(d)-b )d, 

AMN) ~ ^ k\(N-k)\ ' (3 ' 7) 

and the asymptotic value of h(d) should be b\. We apply the same Richardson extrap- 
olation method to h(d) and we compare it with the macroscopic black hole predictions. 
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Figure 1: Microscopic data for f(d) (□), and the Richardson transforms A(d, 2) (a), 
A(d, 3) (o), and A(d, 4) (*). The straight line corresponds to the macroscopic prediction 
b = For the quintic this value is b ps 1.359 and for the available degree 14 the 

Richardson transforms lie 1.8, 2.1, 1.2 % from the macroscopic prediction. For the bi- 
cubic b 0.967, the available degree is higher, 18, and the microscopic counting is within 
.9, 1.2, .3 % from the macroscopic prediction. As an example we give BPS numbers used 
for the analysis at degree 18 of the bi-cubic in table (A.l). 



Calabi-Yau 


dmax 


^4(^maa; 3, 3) 




error 


* 5 (1 5 ) 


14 


1.35306 


1.32461 


2.15 % 


X 6 (l 4 ,2) 


10 


1.75559 


1.71007 


2.66 % 




7 


2.11454 


2.0944 


0.96 % 


A 10 (l 3 ,2,5) 


5 


2.99211 


2.96192 


1.02 % 


^3,s(l 6 ) 


17 


1.00204 


0.987307 


1.49 % 


^4, 2 (1 6 ) 


15 


1.07031 


1.0472 


2.21 % 


^3,2, 2 (1 7 ) 


10 


0.821169 


0.855033 


-3.96 % 


-^2,2,2,2(1 ) 


13 


0.722466 


0.74048 


-2.43 % 


X 4 , 3 (l 5 ,2) 


11 


1.21626 


1.2092 


0.58 % 


X 6 , 2 (l 5 ,3) 


11 


1.52785 


1.48096 


3.17 % 


^4,4(1 4 ,2 2 ) 


7 


1.42401 


1.48096 


-3.85 % 


A6,4(l 3 , 2 2 , 3) 


5 


2.06899 


2.0944 


-1.21 % 




4 


2.95082 


2.96192 


-0.37 % 



Table 3.2: Comparing the extrapolated value of b with the macroscopic prediction. 
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X(6,2) 



A(d,N) 

llr 
10 



Bi-Cubic 



10 



12 
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Figure 2: Microscopic data for /(d) (□), and the Richardson transforms A(d,4) (a), 
A(d,5) (o), and 6) (*). The straight line corresponds to the macroscopic prediction 
b\ = 4 ^|^ . For the degree Xq^ complete intesection this value is b\ ~ 14.44 and for 
the available degree 12 the Richardson transforms lie —11.7, —10.4, —9.77 % below the 
macroscopic prediction. For the bi-cubic bi xs 9.994, the available degree is 18 and the 
microscopic counting is —7.15, —6.88, —6.63 % below the macroscopic prediction. 



Calabi-Yau 


dmax 


A 1 (d max - 3,3) 




error 


estimated 62 


* 5 (1 5 ) 


14 


11.2668 


12.4182 


-9.27 % 


-11.9503 


X 6 (l 4 ,2) 


10 


11.9237 


13.4668 


-11.5 % 


-12.1848 


X 8 (l 4 ,4) 


7 


14.0537 


17.2788 


-18.7 % 


-14.9973 


X 10 (l 3 ,2,5) 


5 


15.2509 


18.8823 


-19.2 % 


-14.9817 


^3, 3 (1 6 ) 


17 


9.29062 


9.99649 


-7.06 % 


-9.63958 


^4, 2 (1 B ) 


15 


10.0226 


10.9956 


-8.85 % 


-10.7834 


-^3,2,2(1 ) 


10 


8.45163 


9.61912 


-12.1 % 


-9.3828 


^2,2,2,2(1 8 ) 


13 


7.84595 


8.88577 


-11.7% 


-8.88773 


X 4 , 3 (T\2) 


11 


9.5981 


10.8828 


-11.8 % 


-9.96404 


X 6 , 2 (l 5 ,3) 


11 


12.5614 


14.4394 


-13.0 % 


-14.2582 


X 4 , 4 (l 4 ,2^) 


7 


9.70091 


11.1072 


-12.7 % 


-9.41295 


X 6 , 4 (l 3 ,2^,3) 


5 


11.1008 


12.5664 


-11.7% 


-10.0821 




4 


11.1378 


12.2179 


-8.84 % 


-8.15739 



Table 3.3: Comparing the extrapolated value with the macroscopic prediction of b\. 
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Two typical examples for the behaviour of the Richardson transforms are plotted in Fig. 
[2j The results for all models are summarized in Table 13.31 

Despite our rather successful verifications of the numerical coefficients bo and b\, we 
should note that the expansion in inverse powers of the charge ( 13.51) is actually an asymp- 
totic series. The asymptotic character of the large charge expansions of microscopic 
degeneracies is manifest in the explicit computations for small black holes in [9] and also 
in the examples we will discuss in section 5. In our case, we can relate the asymptotic 
expansion of (13. 5p to a large genus behavior in a string series, since the coefficients in 
( 13. 5p are proportional to the constant map contribution 

b 9 ~ d g , (3.8) 

where d g is given in (12.111) . This coefficient grows at large g as 

^~(27r)- 4 *(-l)2(2s)!, (3.9) 

which is the typical behavior found in string perturbation theory [SB]- It then follows that 
the series expansion (13.51) for f(d) has zero radius of convergence for any value of d and 
it is rather an asymptotic expansion. Indeed, the d g are the coefficients of the asymptotic 
expansion of the MacMahon function (see [9], Appendix E, for a detailed derivation). For 
these kinds of expansions, the best approximation to their true value (which in this case 
is the function /(d) computed from topological strings) is obtained by truncating the sum 
at the order M which minimizes the error. For an asymptotic series of the form 

oo 

f(w) = W\ h ~ A- k ((3k)\ (3.10) 
fc=i 

the optimal truncation occurs generically at 

JV~4(VVT- ( 3 - U ) 



P\\w\ 

In our case (3 = 2 and we can estimate Af as follows. According to the connection between 
4D/5D black holes [20], the attractor value for the topological string coupling constant is 
g s = Arc |23j. This should be roughly the numerical constant that relates the graviphoton 
field strength to the charge Q in (I2.13p . and it contributes to the coefficients b g an extra 
factor gl 9 ~ 2 ', so that we can refine (13. 8p to 

b g ~ (4vr) 2 %, (3.12) 

and the constant in (I3.10p is A = tt 2 . Therefore, the optimal truncation is at 

A/"~^di (3.13) 
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For the small values of d that we are considering we should therefore expect an optimal 
truncation around M ~ 5 — 10. 

These considerations have implications for our numerical analysis. The Richardson 
method (13.61) is designed in principle for convergent expansions. For asymptotic expan- 
sions, we should expect it to give increasing precision and convergence to the true coeffi- 
cients as long as the order of the transformation N in (13.61) is lower than the truncation 
order M . This is the underlying reason that prevents us from improving the precision 
of the leading coefficients by simply increasing the truncation order N in the Richardson 
method, and we indeed find an oscillating behavior around the expected true value for 
the Richardson transforms with N > 5. 

We try to go one step further and give a rough estimation of the coefficient i>2 in 
(13. ip . which has not been studied in the literature from the supergravity point of view. 
It turns out that the naive method we use for computing the sub-leading coefficient b\ 
gives too big an estimate, which might be a result that the optimal truncation scheme is 
no longer a good approximation at this order. In order to improve this, we use the Pade 
approximation which is well-known for summing divergent series. Given an asymptotic 
series 

oo 

/(*) =$>**, (3-14) 

i=0 

one can evaluate the asymptotic value by defining the following Pade approximation 

p£{z) = ^ i= \? iZ (3.15) 



where the coefficients and Bi are fixed by Taylor expanding the above equation (13.151) 
around z = and match to the first M + N + 1 terms of the original series (I3.14p . 

We take the theoretical values of b and b\ from (13.21) . and use the Monte Carlo 
method to randomly generate the sub-leading coefficients b 2 , ^3 etc, then use the Pade 
approximation to evaluate the asymptotic series (I3.14p for z = ^, where d = 1, 2, • • • , d max . 
We pick the sub-leading coefficients bi (i > 2) that minimize the difference of the Pade 
evaluation with the expected value f(d) from topological strings, i.e. we minimize 



E 



d __, iw - 1)2 (3 - 16) 

We find different values of N, M in the Pade approximation give qualitatively similar 
results. In the last column in Table 13.31 we give the estimated values of sub-sub-leading 
coefficient 62 using the scheme N = 2, M = 1. 

Assuming the constant map contribution is the most significant contribution at this 
order in Q, the coefficient 62 should behave like 

b 2 ~ xk*. (3.17) 
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Figure 3: The plot of —62 vs. 



for 13 Calabi-Yau models. 



We can verify the relation (13.171) by plotting 6 2 against the Euler number x K ® f° r the 13 
Calabi-Yau models we studied. We find as the best fit coefficient 

&2 = 0.047xk*, (3.18) 

see the plot in Figure EJ which is reasonably consistent with the expectation (13.171) . 

From the second row of the Table |3TB"| we can find the numerical values of the genus two 
constant map contribution 62 ~ 0.00017 g gX K ^ ■ Taking into account that g s ~ O(10), this 
is the same order of magnitude as our estimate value of 0.047 from microscopic topological 
string computation. 



3.3 Spinning black holes 



We can try to extract the spin dependence of the black hole entropy from (12.11) . Assuming 
Q ^> J, and expanding in J / Q, we find the following macroscopic prediction for the 
topological string data, 



g m (d) = -^To. 



where 



n(d, 0) 

Q(d, m) 



Po 



3,1- 



(3.19) 



(3.20) 
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d 


(h(d) 


Ai(d, 1) 




A 2 (d, I) 




A*(d, I) 


1 


0.693147 


3.22789 


NA 


NA 


NA 


NA 


2 


1.96052 


6.85432 


NA 


NA 


NA 


NA 


3 


3.59178 


10.9389 


9.03347 


12.2117 


NA 


NA 


4 


5.42856 


14.4696 


9.82804 


13.0403 


12.1257 


6.55334 


5 


7.23677 


16.4156 


10.4705 


12.8183 


11.0112 


10.2148 


6 


8.76658 


16.1819 


10.8618 


11.6135 


10.8785 


9.98996 


7 


9.82591 


13.9173 


10.9692 


9.71239 


10.7516 


8.81357 


8 


10.3373 


10.4832 


10.8121 


7.53259 


10.5093 


7.27017 


9 


10.3535 


7.02869 


10.4477 


5.51774 


10.1494 


5.73628 


10 


10.021 


4.41912 


9.9547 


3.9872 


9.70809 


4.46946 


11 


9.51178 


2.9195 


9.4122 


3.04128 


9.23185 


3.58335 


12 


8.96242 


NA 


8.88129 


NA 


8.76114 


NA 



Table 3.4: The Richardson method for the quintic with spin m — 1, 2, 3. 



For a fixed value m we use again the Richardson extrapolation method to find the 
asymptotic value of g m (d) for large d. We list the values of g m (d) and its first Richardson 
extrapolation A m (d, 1) for spin m = 1, 2, 3, using the quintic as an example. 

We note that the contribution to entropy from angular momentum is proportional to 
d~ 3 / 2 , as compared to the leading static contribution (13. II) of order d 3 ^ 2 . Although the 
prediction (13. 19ft should be the leading spinning contribution, there could be some small 
statistical fluctuation of topological string data which is random for the different spins, 
and which might become comparable to the spinning contribution in (13 . 19[) and result in 
the deviation for large degree d. This can be seen in the quintic example in Table 13.41 
We find that the Richardson series does not converge to an asymptotic value, instead the 
series approach a maximal value before deviating again for large degree d. In order to 
minimize the effect of statistical fluctuation of topological string data, we propose to use 
the extremal values in the Richardson series A m {d, 1) to estimate the asymptotic value 
of po. This is indeed a relatively good estimate for the quintic case where p = 14.9019. 
Other Calabi-Yau manifolds are analyzed in [43] . 

We analyze the 13 Calabi-Yau models using the above approach. Let us define the 
extremal value of the first Richardson extrapolation over the degree d as 

9m = A m {d, 1)| max; (3.21) 

For various Calabi-Yau models and spin m = 1, 2, 3, we compare the value of g m with the 
expected coefficient po given in (I3.20p . The results are summarized in Table [331 We see 
that for larger angular momentum m the deviations become bigger, as expected. 
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Calabi-Yau 


Pn = 3tt(£)5 


Q^ 


gi error 


09 


g2 error 




error 




14.9019 


16.4156 


10.2% 


13.0403 


-12.5% 


10.2148 


-31.5% 


X fi (l 4 ,2) 


11.5429 


12.1492 


5.25% 


10.1828 


-11.8% 


8.21085 


-28.9% 




9.42478 


10.4854 


11.3% 


8.1382 


-13.7% 


5.3473 


-43.3% 


Xmfl 3 , 2,5) 


6.66432 


6.77436 


1.65% 


5.89201 


-11.6% 


3.62439 


-45.6% 


X^(l 6 ) 


19.993 


22.1786 


10.9% 


17.7804 


-11.1% 


14.8114 


-25.9% 


X 42 (l 6 ) 


18.8496 


21.0741 


11.8% 


16.569 


-12.1% 


12.9935 


-31.1% 


9 9 f 1 1 


23.0859 


25.9065 


12.2% 


20.4996 


-11.2% 


16.5636 


-28.3% 


^2,2,2,2(1 8 ) 


26.6573 


30.1999 


13.3% 


23.6923 


-11.1% 


19.2311 


-27.9% 


X 4 , 3 (l 5 ,2) 


16.3242 


17.7685 


8.85% 


14.4772 


-11.3% 


12.2514 


-24.9% 


X 6 , 2 (l 5 ,3) 


13.3286 


15.2332 


14.3% 


11.2819 


-15.4% 


8.06844 


-39.5% 


X 4 , 4 (l 4 ,2 2 ) 


13.3286 


13.9081 


4.35% 


11.618 


-12.8% 


10.6901 


-19.8% 


X 6 , 4 (l 3 ,2 2 ,3) 


9.42478 


9.02611 


-4.23% 


7.87731 


-16.4% 


7.56862 


-19.7% 


X 6>6 (1 2 ,2 2 ,3 2 ) 


6.66432 


5.42333 


-18.6% 


4.91355 


-26.3% 


4.5984 


-31.0% 



Table 3.5: The Richardson method for the 13 Calabi-Yau models with spin m — 1, 2, 3. 



4 Asymptotics of the Donaldson— Thomas invariants 

As we already mentioned, the total free energy of the topological string ( 12.191) can be 
reorganized in terms of Gopakumar-Vafa invariants as in (12.201) . A remarkable property 
of (I2.20p is that for a given class Q G H 2 (X, Z), the expression is exact in the string 
coupling. This is because Castelnuovo's theorem for the ambient space yields n 9 d = for 
d > oty/g for certain a. 

For example, for the quintic the maximal genus g max such that n 9 Q^ ^ fulfills a 
bound 

<W < ^(10 + 5d + d 2 ) (4.1) 

with a decreasing relative deviation in the large d limit. The bound is saturated for curves 
of total degree 5m which are complete intersections of degree (l,5,m) in P 4 , which are 
smooth curves in the quintic. For 5 > m > 1 we can describe the moduli space of the 
D2 brane as follows. The linear constraint has as a parameter space P 4 and allows to 
eliminate one variable from the degree m constraint, which has as many homogeneous 
parameters as monomials in four variables, i.e. as many as there are integer solutions to 
Yli=i n i = m namely ( mH ^ _1 )- The moduli spaces of the curves are therefore fibrations 
of Pv "i ) over P . Using the results of [28J we get for the GV invariant 

n£T = (-l)( m m . (4.2) 

If the bound (14.11) is not saturated for small d the relative deviation can become somewhat 
larger as seen in the Figure HI 
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degree 




Figure 4: Castelnuovo's bound for higher genus curves on the quintic. The dots represent 
n 9max an( j foe curve is (14.11) . 



Let us denote by F'(X,t) the total free energy without the contribution (12.10j) . After 
exponentiation one finds [30] 



Z' GV (X,q,t) = H 
where 



d=l 



oo 2g-2 



g=l 1=0 



>r=l 



q = e 



iA 



(4.3) 



(4.4) 



and we have assumed that there is only one Kahler parameter, so that Q is labeled by a 
single integer d. On the other hand, the conjecture of [35] relating the Donaldson-Thomas 
invariants Dd, n to Gromov-Witten invariants leads to 



Z DT (q,t) = Y J D d,nq n e- dt = Z' Gy (-q,t)M{-q) 



x(x) 



d.n 



where 



X 

n=l v ' 



(4.5) 



(4.6) 



is the MacMahon function. This term reinstalls the constant map contribution. We list 
for reference a few Donaldson-Thomas invariants on the quintic in Table 14. 1 1 

After an extensive discussion of possible tests of the OSV conjecture [37J, the authors 
of [12] isolate as a crucial question for the validity of the latter the growth behaviour of 
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d/n 


-3 


-2 


-1 





1 


2 

















2875 


569250 


1 














609250 


124762875 


2 











609250 


439056375 


76438831000 


3 





8625 


2294250 


4004590375 


1010473893000 


123236265797125 



Table 4.1: Donaldson-Thomas invariants. 



the Donaldson-Thomas invariants. This behaviour is encoded in the scaling exponent /c, 
defined as 

log(£> A 2 d)A3n ) ~ X k . (4.7) 

The question is relevant in the range d 3 — n 2 > for which stable black hole configurations 
exist. 

Because of Castelnouvo's bound, and since our data are up to genus 31, we can calcu- 
late the Donaldson-Thomas invariants exactly in the range < d < 15 and for arbitrary 
high n for the quintic. We are interested in the limit 

k = hm lQglQ ^f^ 3 -' . (4.8) 

A^oo log A 

In order to evaluate it for given values (d, n) we chose A so that d + I = X 2 d for d, I £ N 
and use the fact that log \ D d+ ^ n \ for fixed d, k scales in good approximation linearly with 

n to calculate the interpolated value of the D^+iy at n' = X 3 (d, l)n, with X(d, I) = \J~^-- 
For (d, 0) the latter interpolation is of course completely irrelevant and for charges for 
which the n! values become large it is not very relevant. 

The leading correction to (14.81) is of order l/(log(A)). It makes therefore sense to elim- 
inate this leading correction by logarithmic Richardson-Thomas transforms. We define 

, (0) _ log log |-PA(rf,Q2rf,A(rf,Q3n| Q , 

Kl ~ log X(d,l) ' 1 j 

and the mth logarithmic Richardson-Thomas transform as 

(m) _ k^hogji+^-kt-^iogji) 

1 ~ Mi + 1) - mo • (4 - iUj 

With our knowledge of the topological string up to g = 31 for the quintic we can evaluate 
the Donaldson-Thomas invariants up to degree 15. We plot in the first two graphs <^ 
the data for the kf and its first two logarithmic Richardson-Thomas transforms. The 
graphs clearly indicate that the convergence is improved by the transform. So even if 
there are subleading terms of other forms, we certainly managed to supress the leading 
correction and speed up the convergence. The data further show that there is an universal 
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(1,0) 



2.25 2.5 2.75 3 3.25 3.5 



A 



7 - 
6 
5 
4 

3 - 
2 - 

1 . 4 



(2,0) 



1.6 



2.2 



2.4 



A 



Figure 5: Scaling data k^ (□) and the transforms k^ (a), k^ (o) for the Donaldson- 
Thomas invariants on the quintic in P 4 starting for (d, 0) states. 



behaviour independent of d and that the value of k is within the 2% range close to 2. The 
higher logarithmic Richardson- Thomas transforms are consistent with this value but do 
not determine it better as we also have to take into account values with smaller I hence 
smaller A. We next test the universality of these results for other charges (d, n) in Fig. [61 If 
n/Owe need the interpolation for the n' values. This introduces some random subleading 
errors, which are of the order of the improvement by the second logarithmic Richardson- 
Thomas transform. However as in the figure for (2, 0) we see that higher d seems to lower 
the coefficient of the sub-sub-leading correction and makes already the second Richardson 
transform to converge reasonably well -well enough at least to conclude that the k is 
considerably lower then 3 and very well compatible with the value k = 2 ± 0.03 found for 
the previous charges. We solved the bicubic in P 5 up to genus 29, which yields complete 
informations about the Donaldson-Thomas invariants up to degree 18. A similar analysis 
as above confirms the analysis for the quintic. The corresponding plots are in Fig. [7] and 
Fig. [HJ Again a detailed summary of the data for more models can be found at |43j . 
We note a slight noise in the transform k^ in Fig. [HI which is presumably due to the 
interpolation in the n value of D^^ described above. The results for the other models are 
similar, but somewhat less precise due to smaller values of d that are currently available. 

To summarize: our analysis indicates that the value of k is indeed universal and close 
to k = 2. This strongly suggests that the "mysterious cancellations" [12] that eventually 
make possible to extend the the OSV conjecture to small coupling, actually take place. 
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1 . 4 



1.6 



2.2 



A 



2 - 
1.2 
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1.! 



-A 



Figure 6: Scaling data (□) and the first transform A;* 1 ) (a) for the Donaldson-Thomas 
invariants on the quintic in P 4 for the (2, 1) and (3, 1) states. 
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Figure 7: Scaling data k^ (□) and the transforms A;' 1 ) (a), k^ (o) for the Donaldson- 
Thomas invariants on the bic-cubic complete intersection in P 5 starting for (d, 0) states. 
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Figure 8: Scaling data £j(°) (□) and the first transform A;* 1 ) (a) for the Donaldson-Thomas 
invariants on the bic-cubic complete intersection in P 5 for the (3, 1) and (3, 2) states. 
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5 K3 fibrations 



5.1 Topological strings on K3 fibrations 

We will now consider Calabi-Yau manifolds X that have the structure of a K3 fibration, 
i.e. there is a fibration of the form 



7T 



:X->P\ (5.1) 



where the fibers are K3 surfaces. When the fibration is regular the homology of X can 
be written as 

H 2 {X, Z) = ([P 1 ]) © Pic(K3), (5.2) 

where Pic(K3) is the Picard lattice of the K3 fiber. The rank of this lattice will be denoted 
by p, and £ a , a = 1, • • ■ , p will denote a basis for this lattice. Let uj be the complexified 
Kahler form on X. The complexified Kahler parameters of X are given by 

S = I co, t a = I u, a=l,---,p. (5.3) 

We will denote by r]s, i] a the two-forms which are dual to P 1 , S a . 

It turns out that type IIA string theory compactified on these manifolds is very often 
dual to heterotic string theory compactified on K3xT 2 [271 ED]- Under this duality, 5* 
becomes the axidilaton of the heterotic string. It follows that in the regime S — > oo one can 
map computations in the type IIA theory to perturbative computations in the heterotic 
string. In particular, the F g couplings of topological string theory (which are graviphoton 
couplings in type IIA theory) can be computed exactly at one-loop in the heterotic string, 
provided the Kahler parameters are restricted to the K3 fiber J27J [311 [251 12] • We will now 
review here some of these results. 

The topological string amplitudes F g (S, t) on these fibrations have the following struc- 
ture, 

F (S, t) = \c abc t a t h t c + \c ab St a t h + ^x(X) + ^ (t) + 0(e~ s ), 
biz 

Fi{S, t) = ^(c s S + cX) + + 0(e~ s ), ^ 

F g (S, t) = d gX (X) + F g (t) + 0(e~ s ), g>2. 

In these formulae, C a bc and C ab are triple intersection numbers in the fiber and in the 
mixed fiber/base direction, respectively. Notice that 

Cabc = / Vahr] b A r) c , C ab = / r)s A i] a A (5.5) 
Jx Jx 

We also have 

c a = J c 2 (X)Ar] a , a = 1, ■ ■ ■ , p, c s = J c 2 (X) A r} S . (5.6) 



22 



For K3 fibrations with trivial fundamental group one has cs = 24 [36] , but for the Enriques 
Calabi-Yau (which we will also analyze), eg = 12. The coefficient d g is the contribution 
of constant maps written down in (12. lip . In (j5.4j) . T g (f) denotes the contribution of 
worldsheet instantons in the K3 fiber. It follows from [21 ESI ESI EH ED] that the F g {t) can 
be completely determined in terms of a single modular form that we will denote fx(l)- in 
order to write down an explicit formula for J- g (t), we have to introduce the quasimodular 
forms V g (q) which are defined by 

(S^J-t^r'VM (5.7) 

The quantities V g (q) can be explicitly written in terms of generalized Eisenstein series 
[34] . and one has for example 

Vi(q) = ^E 2 (q), V 2 (q) = ^( 5E i + (5.8) 
We now introduce the coefficients c g (n) through 

V g (q)fx(q) = J2 c f^ n - ( 5 - 9 ) 

n 

One then has the following expression for the heterotic J~ g {t): 

= E cf(g 2 /2)Li 3 - 2;? (e-^), (5.10) 

QePic(K3) 

where Li n is the polylogarithm of index n 

OO j, 

fc=i 

In (15.101) we have also denoted 

Q-t = n a t a } Q 2 = C ab n a n b} (5.12) 

where C ab = is the intersection form of the Picard lattice Pic(K3). 

We will particularly interested in three special K3 fibrations: the STU model, the ST 
model, and the Enriques Calabi-Yau. Let us give some extra details for these cases: 

• The STU model has p = 2 and it can be realized by a complete intersection in a 
weighted projective space which is frequently denoted by X 2 4(l, 1,2,8, 12). It has 
Euler characteristic x = — 480. The classical prepotential can be obtained from the 
nonvanishing intersection numbers, 

Cm = 8, Cn2 = 2, Cn = 2, C12 = 1, (5.13) 
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while the classical part of Fi(S,t) is encoded by 

ci = 92, c 2 = c s = 24. (5.14) 

The modular form encoding the information about topological string amplitudes in 
the fiber is given by [34J 

/s T u(g) = -^5T^)- (5-15) 
It is sometimes useful to parametrize the Kahler cone in terms of the variables 

T = t 1 +t 2 , U = tx, (5.16) 
In this basis one has Q 2 /2 = mn. 

The ST model has p = 1 and is realized in type IIA by the CY Xu(l, 1, 2, 2, 6). It 
has x = — 252 and the classical intersection numbers 

Cm = 4, C n = 2, (5.17) 

as well as 

Cl = 52, c s = 24. (5.18) 
The Kahler parameter along the fiber is usually denoted as 

T = t 1 . (5.19) 

The relevant modular form is 1291 [301 



Arfa) = -^(9), (5.20) 



where 

%) = E^=^(r/2), 

F> ^ ^J(r/2). (5 ' 21) 
F 6 = E 6 - 2F 2 (9 4 - 2F 2 )(6 A - 16F 2 ). 
Notice that Q 2 /2 = ra 2 /4. 

The Enriques Calabi-Yau is given by the free quotient (K3 x T 2 )/Z 2 , and was 
introduced in the context of type II/heterotic duality in [18]. It is an elliptic fibration 
with p = 10. It has C a b c = 0, while C a b is given by the intersection numbers of the 
Enriques surface E. The Picard lattice is 

Pic(K3) =r 1 ' 1 ©E 8 (-l), (5.22) 
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and 

c a = 0, c s = 12. (5.23) 

The topological string amplitudes in the fiber were obtained in [32j (see also [22]). 
They are also controlled by a single modular form 

Uq) = ~^wy (5 - 24) 

but their form is slightly different from (I5.10p 

F g {t)= E c 9 £ (Q 2 )(2 3 -^Li3_ 2s (e-«-VLi3- 2ff (e-^; 

Q£Pic(K3) ^ 

where c^(n) are defined again by (I5.9p . 



(5.25) 



5.2 Microscopic degeneracies and their asymptotic expansion 

We have seen that, at least in the case of topological strings on K3 fibrations, and for 
classes Q restricted to the K3 fiber, one can obtain closed formula for the topological 
string amplitudes at all genera. It should be therefore possible to extract a closed formula 
for the generating functional of Gopakumar-Vafa invariants. In fact, by using the product 
formula 

oo 

■& x {v\t) = -2g^sin(7rz/) JJ(l-g")(l-2cos(27rz/)g" + g 2n ) (5.26) 

71=1 

one finds from the expression (15.101) and the structure (12.201) 

oo 

E E^ 2 > Q2/2 = fx(p)e( Z ), (5.27) 

QePic(K3) r=0 

where l;(z) is the function that appears in helicity supertraces, 

M = II (l-p*)2 + Z 2pn = II {1 - p n y){1 _ p n y -iy ^ 

where we have set z = —i(y^ — y~^)- 

We can now obtain a closed formula for the microscopic degeneracies. In order to have 
a description which incorporates as well the elliptic genus, we will count the microstates 
as in (12.181) but with r — > r — 1. With this definition, the l.h.s. of ( 15.271) . expanded in q, y, 
is precisely the generating function of microscopic degeneracies Q(Q, m), summed over all 
m,Q. We then arrive to the expression 



oo 

E E 

QSPic(K3) m=-oo 



tt{Q,m)y m p Q2/2 = /x(p)e>,*), (5-29) 
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where we have written 

y = e 2 ™ p = e 2 ™. (5.30) 

Notice that if we consider X = K3 x T 2 and restrict to classes Q in the fiber, the counting 
of microstates given by the elliptic genus is 

OO j 

X (S P K3; q, y) q0 = ]J {1 _ pNyo{1 _ p N y) 2 {1 _ pNy -iy = ^ff^ ^ 



N=l 



This has the same form than ( 15.29R with 

/ K3 xT 2 (p) = ^y, (5.32) 

therefore we can consider the "small" D1-D5 system as a particular case of our analysis. 

The expression (15.291) tells us that the microscopic degeneracies we are looking for are 
simply the Fourier coefficients of the object in the r.hs. We can then invert it to write 

/•|+i0+ pi 

n(N,m)= da due- 2 ^ Na+mu) ^(u } a) } N = Q 2 /2, (5.33) 

</-i+i0+ Jo 



where we defined 

$(v,a) = f x (p)e(iy,a). (5.34) 

and we have assumed that N is a non-negative integer (this can be guaranteed by rescaling 
p — > p k for some appropriate k). The contour in (15.331) has been chosen to avoid the poles 
in the integrand. 

We will now evalute the asymptotic expansion of Q(N) = Q(N, 0) in inverse powers 
of N. Nonzero values of the spin m = can be analyzed in a similar way. The expression 
we will find is exact up to corrections which are exponentially suppressed in the large 
charge limit iV — >• oo. Notice that in our situation we can not appeal to the Rademacher 
expansion which was used in [T5J [H] , since (I5.34p is not a Jacobi form (it can be regarded 
as a Jacobi form with negative index). It is likely that an analog of the Rademacher 
expansion exists, but we will perform a direct evaluation of the integral f!5.33j) in the 
spirit of the counting of states with spin in Appendix C of j9] and in [TO] . 

First of all, we reexpress the integrand (15.341) in terms of di{v\a) as, 

${u,a) =4sin 2 (vr^ 6 (p)-^M Y . (5.35) 

Using the modular behavior of $i(i/\a) under the S transformation a — > a = — 1/er we 
get, 

M ^ a) = -J^ e ^H) sin (™)\i + o(e-^)). (5.36) 
V— icr V o J I J 
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It is easy to see that the saddle point evaluation of (15.331) is governed by 

'• = ^ +0 Q- (5 - 37) 

Therefore, the corrections to (1 5 . 3 6 [) will be exponentially suppressed. Using the modularity 
of rj{p), and taking the part of the sin in the denominator which is not exponentially 
suppressed, we obtain, 

$(z/, a) ~ -4a- 2 e 2 ^ u2 ~ u) sin 2 (vri/)/ x (p). (5.38) 
Therefore, in order to compute the asymptotics of ( 15.331) we just need 



fi(iV)~-4/'" 1 dae- 2 ^^- f dve 2 ^ 2 -^sm 2 (7rv). (5.39) 

J-l+iO+ O" Jo 

The integral over v is easily worked out in terms of the error function Erf (a;), as follows, 

+ yi^ , { Erf (& +i, )- Erf (&- i) )}- 

Due to (15.371) we can use the asymptotic expansion of the Erf function, 

Erf(x) ~ 1 - V(-l) r l J " x- (2r+1) , \x\ -> oo, |arg(-x)| < tt. (5.41 



(5.40) 



7T 

r=0 



Ignoring terms which are exponentially suppressed at large N, we find, 

/■l 1 00 ;l+3r 

/ dz/e 2 ^ 2 -^sin 2 (7rz/)~--^^(2r-l)!!G r (^), (5.42) 

^0 n ^ 



r=0 



with, 

Again, due to (I5.37p . we can expand it around a = 0, 



^)^E( 2(i+ r r) )' 

„_ n \ / 



(7 3+2s+r . (5.44) 



s=0 

Putting all together, we obtain, 

\+i0+ 



r=0 s=0 v '2" 
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We now work out the integral, 

/•|+i0+ 

A sr {N)= dae- 27rlN<7 f x {p)a 1+2s+r . (5.46) 

We assume that fx(p) has modular weight w, so that fxip) — cr ~ w fx(p), where p = e~~ . 
For the modular forms that we consider here, fx{P) — cp~ a + ■ ■ • , and the integral above 
gives a modified Bessel function 

A s , r (N) ~ ci 1+2s+ ^4 +r+2 _ w (47rV^V). (5.47) 

We end up then with the following result for the exact asymptotics of the microscopic 
black hole degeneracy, 

n(JV) ~ 2ci- f; f2(-iy ( 2(1 + * + r) ) W,^(W^V). (5.48) 

r=0 s=0 ^ ' 

Using now the formula for the asymptotic expansion of I functions (see for example App. 
A of [9]), we find for the entropy S(N) = log Q(N) the following expansion 

S ~ 47TV/OJV — log iV + log + — -= = + 0(N ). 

4 V 7r / 32n^a y/N 



(5.49) 

The expansion in powers of l/iVa in (15.481) . which is obtained by using the asymptotics 
of modified Bessel functions, is the expansion of the original integral around the saddle 
point (|5.37p . This can be verified by an explicit computation of the first few orders of the 
saddlepoint expansion. 

Let us now evaluate (I5.49p in some examples. For K3xT 2 we have (w,a,c) = 
(—12, 1, 1), and the entropy reads 

29, / a r\ i i fV%\ 591 1 , ^ /Ar _!, 



S ~ AnVN - — log(iV) + log! ^- 1 - — -j= + 0(N- L ). (5.50) 



(5.51) 



For the STU model, with the values (w, a, c) = (—2, 1, —2), we find 

, q 19Q 1 

S ~ 4nVN - - log(iV) + log ^ + -1— + 0{N- 1 ). 
4 V 77 / 3271-ViV 



The ST model is slightly different, since in /st(p) both integer and rational powers of 
p appear. As mentioned above, we should redefine p — > p 4 and write down the generating 
functional for the degeneracies as 

oo 

E m,m)y m p 2Q2 = / ST (p 4 )£>,4a), (5.52) 

QePic(K3) m=-oo 
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where we recall that M = 2Q 2 = n 2 is an integer. The asymptotics is given by the integral 
n ST (M)~- dae -^MahA^l /■ 1 dvsin 2( 7rv ) e ^(^-). (5.53) 

J-l+iO+ ^ Jo 

The integral over v is given by (I5.40p upon replacing a — ► 4a. Since, 

W) = -2^^p^ 3 (2 ( r) ~ -IGv^ale^, (5.54) 



one finds in the end, 

, +2s+r (27rVM), (5.55) 



~,_n s— n \ / 



r=0 s=0 

and from here one can read the entropy, 



S(Q) ~ 47r^g 2 - 2 log(Q 2 ) + • • • (5.56) 

Finally we turn to the case of Enriques CY manifold. It follows from (15.251) that one 
has to distinguish two types of homology classes: the classes Q whose entries contain at 
least an odd integer (which were called odd classes in [32]), and the classes Q for which 
all entries are even (called even classes). A simple calculation shows that the generating 
function of Gopakumar-Vafa invariants for the odd classes is given by 



E E ««P° 2 ^ 1 = i^7k(« 2 (''Art-«'V2.-rt)- 

r=0 Q odd * bm I 2 ) 



f E {q) 

r=0 Q odd 

while for the even classes is given by 



(5.57) 



S £ n^^^ = --^) (e 2 W2,p)-f(,/2 ) -p)) 

r=0Qeven 4 Sm \T ) (5.58) 



Notice that for even classes Q 2 = mod 4, while for odd classes one only has Q 2 = 
mod 2. In contrast to the previous K3 fibrations, in the above generating function we 
have , instead of p® / 2 , and this will lead to a different leading term as compared for 
example to the STU model. 

The computation of the asymptotics of the microstates is similar to the one that we 
just performed. Let us begin with odd classes. Using the identity, 



^=^ %iXSly (5 ' 59) 
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and proceeding as in the previous case, we find, 

n odd (N) = n^N) + n 2 (N), n = q 2 /2, (5.60) 

where, 

4+io+ 



dae-^ iN(T a 2 r ] 6 {2a)^l{2a)f E (p) / dz/sin 2 (vrz/)e^^ 2 -^ + ^, 

-r° + + ' ° (5.61) 

Q 2 (N)~-4i dae- 47TiNa aT] 6 {a)f E {p) dz/sin 2 (7rz/)e^ (i/ - 1)2 . 

J-~+iO+ JO 

As before, we evaluate the integrals over v in terms of the Erf function and its asymptotic 
expansion. We then use the modularity properties of the different functions involved here 
to obtain, 

1 " (2r-l)!! ^ /2(l + a + r)V ,_ (1+r+s) , f 

i — n o— n V / 



(5.62) 

Let us now consider the even classes, (I5.58p . Comparing (15.581) with (I5.57p . we see that, 

n even (N) = n odd (N) - h(N) (5.63) 

where, 

/I+iO+ „l 
da dz/e- 4i7rAr(J 4sin 2 (7rz/)/ s (p 4 )(e 2 (^p 4 )-e 2 (^ V))- (5.64) 
i+iO+ Jo v J 

A computation similar to the one we performed shows that Q(N) is exponentially sup- 
pressed with respect to Q odd (N), since it leads to terms that go like exp(-K\ / 2N) and 
exp(ir\/QN). Therefore, as an asymptotic expansion in 1/y/N, £l even (N) ~ f2 odd (iV), and 
the asymptotics does not distinguish between the even and the odd classes. We finally 
obtain, for the small Enriques black hole, 

S E (Q) ~27r v ^-ylog v / Q 2 + --- • (5.65) 

The main conclusion of our analysis is that, in all cases, the leading term of the 
microscopic entropy for these black holes is given by 



5(g)~27r^Q 2 , (5.66) 

since c$ = 24 for K3xT 2 , the STU and the ST models, but c$ = 12 for the Enriques CY. 
Of course, our analysis has also given precise formulae for the subleading terms. 
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The leading behavior (I5.66P can be also verified by a numerical analysis similar to the 
one performed in sections 3 and 4. For example, for the STU model we have computed 
the quantity f(N) = S(N)/y/N for 1 < iV < 50, where S(N) = log Q(N). In order 
to subtract the logarithmic term in the asymptotic expansion (I5.5ip we consider the 
t i*ctiisf orixi 

A(N) + 1)S(N + 2) - (2N + 1)S(N + 1) + NS(N) 

(N + l)\/W+~2 - (2N + l)i/FTI + Ny/N 

In Fig. [9] we plot f(N) (bottom) and A(N) (top). The horizontal line is the expected 
asymptotic value An for both quantities as N — > oo. As before, the transform A(N) 
improves rapidly the convergence. 
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Figure 9: Microscopic data for f(N) = S(N)/vN (bottom) and its transform A(N) 
(top), defined in (I5.67p . for the STU model, and for 1 < iV < 50. The horizontal line is 
the expected asymptotic value 4n. 



5.3 Macroscopic entropy for small black holes 

The 5d black holes obtained by wrapping the M2 branes along cycles in the K3 fiber have 
actually vanishing classical entropy and are therefore small black holes. Indeed, as we 
have seen, the leading asymptotic degeneracy scales like Q, and not like Q 3 ^ 2 . This is also 
what is found for small 4d black holes [9]. 

Let us briefly show that the classical area of these black holes is zero for any set of 
intersection numbers C a b c , C^. In order to do this, we can use the 5d attractor mechanism 
described in section 2. Equivalently, by using the 4d/5d connection of [20], we can map 
the 5d black hole to a 4d black hole with D6 charge p° — 1 and D2 charges Qa- At the 
level of the leading macroscopic entropy, the 4d computation gives the same result as the 
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5d computation [20J. In the 4d language, we start with the tree level SUGRA prepotential 

1 X s X«X b 1 X*X^ 
- ~2 — Q Gabc — — • 

We will do the computation for a generic D6-D2 charge, i.e. we will start with generic 
charges p°, Q a , Qs, and then take the charge Qs — > at the end of the computation 
(as well as setting p° = 1). This will guarantee that we obtain generic solutions to the 
attractor mechanism. 

Let us first assume that C abc = 0, as it happens in K3xT 2 and the Enriques Calabi- 
Yau. In this case, the attractor equations are easily solved as, 



(X2,X*,X;) = [p^H ^,hM^Q a I (5-09) 



, zQs'l Q 2 

where 

Q 2 = C ab Q a Q b , Q a = C ab Q b . (5.70) 

The entropy is given by 

S = 7r^2p°Q s Q 2 , (5.71) 

and it vanishes in the limit Qs — > 0. This is as expected. 

If we now consider a general prepotential with nonvanishing C abC) the attractor equa- 
tions are now solved at 



(X°, X s , X a ) = (p° , i^ ^-J s , iy/^Q^j , (5.72) 
where the £ A are solutions to, 

£ a f = 1 

« , , (5-73) 

Notice that, in these variables, the model with C abc = corresponds to the smooth values, 

={^Q^-%)- ( 5 - 74 ) 



/Q 2 

We can already see that, in the limit Qs — > 0, the perturbation by C abc in (15.731) vanishes, 
therefore in the limit of zero charge in the base the presence of nontrivial intersection 
numbers in the fiber should be unimportant. More formally, it is easy to see that one can 
construct a consistent solution of (15.731) of the form, 

oo 
n=l 

32 



where the coefficients depend on C a bc and can be calculated order by order. In terms 
of the £ the macroscopic entropy is 

s = n^2fQ~s(c ab eee + \Q s c a u a e^ , (5.76) 

and, in the limit Qs — > 0, it will vanish irrespectively of the value of C a b c . Therefore, 
5d black holes whose membrane charge is restricted to the K3 fiber of a K3 fibration are 
always small. This can be checked as well by detailed computations in different models 
(like the STU and ST models considered above). 

Since the leading contribution to the entropy vanishes we should now look at the 
sub leading terms in the macroscopic entropy. As we explained in section 2, it was shown 
in [231 [U [8] that these terms are obtained by performing the shift 

Qa — Q A = Qa + (c 2A , C = o ■ ( 5 ' 7T ) 

o 

The leading term in the entropy for the small 5d black hole is given (for large charge Q) 
by performing this shift in (15.711) 



S = 2k^Q*. (5.78) 

This can be derived in detail by solving the attractor equations with shifted charges (15.771) 
as a power series in 1/Q, and then taking the limit Qs — > 0. Notice that the entropy 
(15.781) only depends on C a b and cs- Also, in this regime, the solutions of the attractor 
equations occur at values of the Kahler parameters which are of the order of the string 
size, and the SUGRA calculation might be problematic. Indeed, it is easy to see that 
(15.781) does not agree with the leading term of the asymptotics that we obtained in the 
previous subsection. By comparing (I5.66P with (15.781) we find that the formula agree if 
we set instead £ = 1/6. This is the value of ( that is predicted by the 4d/5d connection 
of [20]. 

In [2"3"j [S] it was noticed that the subleading correction (15.771) obtained in a macroscopic 
5d computation was not in accord with the subleading correction predicted by [2U] and 
the 4d attractor mechanism. We now find that, for big 5d black holes, the subleading 
correction for the microscopic entropy is in rough agreement with (15.771) . while for small 
5d black holes the leading asymptotics is in accord with a 4d computation for a small 
D6/D2 system with p° = 1. As we already mentioned, in the case of small black holes, 
the SUGRA computations with which we are comparing our results should receive large 
corrections, but in other situations they still lead to results which are in agreement with 
the microscopic counting, as in [9l [13] . In our case we obtain a result in disagreement with 
the 5d computation but in agreement with the 4d computation. It would be interesting 
to resolve this puzzle. 
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6 Conclusions 



In this paper we have studied the microscopic counting of 5d black hole states by using 
topological string theory. In the case of big black holes, we have given convincing nu- 
merical evidence that the BPS invariants encoded in the topological string amplitudes 
account correctly for the macroscopic entropy of spinning black holes. Moreover, we have 
also shown that the data favour the "mysterious cancellation" of [12J that makes possible 
to extend the validity of the OSV conjecture, and we were able to explore new aspects 
of black hole entropy which have not been studied before using supergravity. Clearly, it 
would be very desirable to improve our numerical results with more data. Using the inter- 
play between modularity and an-holomophicity in topological string theory [I2J [26j [22] . 
analytic results on the asymtotics might be not out of reacffl 

We also gave exact formulae for microscopic degeneracies of a class of small 5d black 
holes, which are obtained by wrapping M2 branes in the fiber of a K3 fibration, and we 
computed the asymptotic expansion in inverse powers of the charge. As expected, the 
calculation shows that for small black holes the leading term in the entropy scales like 
S — ► XS when the charges are scaled with A. We found however that the coefficient of 
the leading term does not agree with the shift of charges obtained in [23], [1] [8] in a 5d 
SUGRA computation. In principle there is no reason why these two computations should 
agree, since small M = 1 black holes are generically beyond the SUGRA approximation. 
On the other hand, the microscopic results are well reproduced by the 4d/5d connection 
of [20] and a 4d attractor computation. We should emphasize however that for big black 
holes the 5d shift (I5.77P fits our data better than the 4d shift with £ = 1/6. It would be 
very interesting to understand this better. 
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A General features of the instanton expansion 

The asymptotic behaviour at the conifold, Castelnouvo's theory, and the calculation via 
degenerate Jacobians, suggest some general features of the Gopakumar-Vafa expansion. 

4 Recently beautiful analytic proofs of the asymptotic of the Fourier coefficents of Mock-Theta functions 
have been obtained using a somewhat similar interplay [6]. 
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Our data for the 13 one-parameter models suggest further universal features. The purpose 
of this appendix is to describe some of these general features. Typical data for high degree 
look as is table IA.1I 

The last nonzero entry is from the smooth genus 28 complete intersection curv^f] 
(1,2,3,3) of degree 18. By Castelnouvo's theory g = 28 is the largest possible genus for 
degree 18. The degree one constraint parametrizes an P 5 . The moduli space AiH is a 
fibration of this P 5 over a projectivization of the 15 parameters in the quadratic constraint. 
I.e. MH is the total s P ace of p5 pM > with Euler number x(M 9 d Zf 8 ) = 5 x 15 = 90 and 
nil = (-1) 5 + 14 90 = -90. 

As it can further be seen in table \A.1\ the numbers grow from genus g = to g = 3 
and fall thereafter. This feature might be related to the binomials in the description of 
the moduli of space as a singular fibration of the Jacobian Jac28 of the g = 28 curve over 
Aifl- In this description the contribution of a g = g — 5 curve comes from degenerating 
the genus 28 curve with 5 nodes. As explained in [28J the contribution of the degenerate 
Jacobians can be expressed by the Euler numbers of relative Hilbert schemes C^ n ' as 

5 

n 9 i -s^ { _ 1)d un(M)+6J2b(g-p,8-p)x(C^) , (A.l) 

with b{g,k) = ( 2g ; 2 ). A simple G auss approximation of binomials fits the behaviour of 
the n d for large d relatively well. We show this in Fig. [10] for the bi-cubic at degree 27. 
The numbers n d are exact and in contrast to (1A.1I) they count correctly all contribution 
from colliding nodes, all contributions from reducible curves as well as contributionsfrom 
smooth curves in the class d with genus g < g. 

Very important for the cancellations in the asymptotic behaviour of the Donaldson- 
Thomas invariants is the occurrence of negative numbers. While it is clear that such 
contributions can arise if the dimensions of the D-brane moduli space is odd, we do 
not understand a priori the remarkable pattern with which these signs occur. The first 
occurrence of negative signs at g neg (d) is graphed for the quintic and the bi-cubic in 
Fig. [TTJ The data suggest that g neg (d) follows a parabola similar to the Castelnouvo 
bound. From the first occurrence of the negative sign the n 9 d are alternating in sign for 
g <C g. For g ~ g the behaviour becomes more erratic. The Gauss approximation for the 
absolute values of the n d and the sign pattern is very characteristic of the degeneracies of 
microstates of a large black hole. In contrast the absolute value of the n d is falling and 
the signs are alternating with (— l) 9 starting at g = for small black holes as shown for 
example for the ST-model. 

5 A complete intersection curve (1, n, 3, 3) with degree 9n has in general genus g = |(1 + 3n)(2 + 3n). 
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genus 


degree=18 





144519433563613558831955702896560953425168536 


1 


491072999366775380563679351560645501635639768 


2 


826174252151264912119312534610591771196950790 


3 


866926806132431852753964702674971915498281822 


4 


615435297199681525899637421881792737142210818 


5 


306990865721034647278623907242165669760227036 


6 


109595627988957833331561270319881002336580306 


7 

8 


28194037369451582477359532618813777554049181 
5218039400008253051676616144507889426439522 


9 


688420182008315508949294448691625391986722 


10 


63643238054805218781380099115461663133366 


11 


4014173958414661941560901089814730394394 


12 


166042973567223836846220100958626775040 


13 


4251016225583560366557404369102516880 


14 


61866623134961248577174813332459314 


15 


451921104578426954609500841974284 


16 


1376282769657332936819380514604 


17 


1 186440856873180536456549027 


18 


2671678502308714457564208 


19 


-59940727111744696730418 


20 


1071660810859451933436 


21 


-13279442359884883893 


22 


101088966935254518, 


23 


-372702765685392 


24 


338860808028 


25 


23305068 


26 


-120186 


27 


-5220 


28 


-90 


29 






Table A.l: Gopakumar-Vafa invariants n g d in the class d — 18 for the complete intersection 
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Figure 10: The binomials dominate the behaviour of large d Gopakumar-Vafa invariants. 
For the degree 27 class on the bi-cubic we find n| 7 ~ e i67.747 e -0.0985(s-9.i08) 2 



9 


d= 1 


2 


3 


4 


5 


6 


7 





2496 


223752 


38637504 


9100224984 


2557481027520 


805628041231176 




1 





-492 


-1465984 


-1042943520 


-595277880960 


-316194812546140 




2 





-6 


7488 


50181180 


72485905344 


70378651228338 




3 











-902328 


-5359699200 


-10869145571844 




4 











1164 


228623232 


1208179411278 




5 











12 


-4527744 


-94913775180 




6 














17472 


4964693862 




7 
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A further remarkable fact is the very universal scaling for the maximal value M(d) for 
n 9 d for given d. This value behaves like 

M{d) = exp ((a + 6d) 4/3 ) (A.2) 

with very similar values for a and b for different one-parameter models, as shown for the 
quintic and the bi-cubic in Fig. [121 
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gneg(d) gneg(d) 




d 



Figure 11: The first occurrence of negative n g d for the quintic (on the right) and the bi- 
cubic (on the left). The fit is m(d) = a + bd + cd 2 with a = —4.6, b = .94 and c = .019 as 
well as a = —5.2, b = 1.0 and c = .017 for these two , respectively. 




Figure 12: m(d) = log(M(<f)) 3 / 4 for the quintic on the right and the bi-cubic on the left. 
a = 5.164 and b = 1.511 as well as a = 5.202 and b = 1.509 for the cases plotted. 
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